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1 Introduction 

A general problem in physics is to deduce the macroscopic properties of a quan- 
tum system from a microscopic description. Such systems can only be described 
mathematically on a scale much smaller than the scales which are probed exper- 
imentally or on which the system naturally interacts with its environment. An 
obvious reason is that systems consist of particles whose individual behaviour 
is known and also whose interactions with neighbouring particles is known. On 
the other hand experimental probes interact only with systems containing large 
numbers of particles and the apparatus only responds to their large scale aver- 
age behaviour. Statistical mechanics was developed expressly to deal with this 
problem but, of course, only provides a framework in which detailed methods of 
calculation and analysis can be evolved. 

These notes are concerned with the physics of phase transitions: the phenomenon 
that in particular environments, quantified by particular values of external param- 
eters such as temperature, magnetic field etc., many systems exhibit singularities 
in the thermodynamic variables which best describe the macroscopic state of the 
system. For example: 

(i) the boiling of a liquid. There is a discontinuity in the entropy. 



where AQ is the latent heat. This is a first order transition; 

(ii) the transition from paramagnetic to ferromagnetic behaviour of iron at 
the Curie temperature. Near the transition the system exhibits large-range 
cooperative behaviour on a scale much larger than the inter-atomic distance. 
This is an example of a second order, or continuous, transition. Scattering 
of radiation by systems at or near such a transition is anomalously large 
and is called critical opalescence. This is because the fiuctuations in the 
atomic positions are correlated on a scale large compared with the spacing 
between neighbouring atoms, and so the radiation scattered by each atom 
is in phase and interferes constructively. 
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Most of the course will be concerned with the analysis of continuous transitions 
but from time to time the nature of first order transitions will be elucidated. 
Continuous transitions come under the heading of critical phenomena. 
Broadly, the discussion will centre on the following area or observations: 

(i) the mathematical relationship between the sets of variables which describe 
the physics of the system on different scales. Each set of variables encodes 
the properties of the system most naturally on the associated scale. If we 
know how to relate different sets then we can deduce the large scale prop- 
erties from the microscopic description. Such mathematical relationships 
are called, loosely, renormalisation group equations, and , even more 
loosely, the relationship of the physics on one scale with that on another is 
described by the renormalisation group. In fact there is no such thing 
as the renormalisation group, but it is really a shorthand for the set of 
ideas which implement the ideas stated above and is best understood in the 
application of these ideas to particular systems. If the description of the 
system is in terms of a field theory then the renormalisation group approach 
includes the idea of the renormalisation of (quantum) field theories and 
the construction of effective field theories; 

(ii) the concept of universality. This is the phenomenon that many systems 
whose microscopic descriptions differ widely nevertheless exhibit the same 
critical behaviour. That is, that near a continuous phase transition the de- 
scriptions of their macroscopic properties coincide in essential details. This 
phenomenon is related to the existence of fixed points of the renormalisation 
group equations. 

(iii) the phenomenon of scaling. The relationship between observables and pa- 
rameters near a phase transition is best described by power-law behaviour. 
Dimensional analysis gives results of this kind but often the dimensions of 
the variables are anomalous. That is, they are different from the obvious 
or "engineering" dimensions. This phenomenon occurs particularly in low 
dimensions and certainly for o? < 4. For example in a ferromagnet at the 
Curie temperature Tc we find 

M ~ h^, 

where M is the magnetization and h is the external magnetic field. Then 
the susceptibility, x = ^! behaves like 

Since 5 > 1, x diverges as ^ — >• 0. The naive prediction for (5 is 3. 5 is 
an example of a critical exponent which must be calculable in a success- 
ful theory. The coefficient of proportionality is the above relations is not 
universal and is not easily calculated. However, in two dimensions the con- 
formal symmetry of the theory at the transition point does allow many 
of these parameters to be calculated as well. We shall not pursue this topic 
in this course. 
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Note. The general ideas of relating physical phenomena on widely differing 
scales by renormalisation group methods is widely applicable in many fields. For 
example: 

(i) general diffusion models such as diffusion on fractal structures and the large 
scale effects of diffusive transport processes in fluid flow; 

(ii) turbulence in fluids. The fluid velocity field, u(x, t), of fully developed 
turbulence has energy density, as a function of wavenumber k. 



with E{\s.) ^ k 3 for large k. 

This result is derived by "naive" dimensional analysis (Kolmogorov). How- 
ever, a full solution of the Navier-Stokes equation can correct the exponent: 



where rj is the anomalous term. 

2 The Phenomenology of Phase Transitions 

Statistical systems in equilibrium are described by macroscopic, thermodynamic, 
observables which are functions of relevant external parameters, e.g., tempera- 
ture, T, pressure, P, magnetic field, h. These parameters are external fields 
(they may be x, t dependent) which influence the system and which are under 
the control of the experimenter, 
the observables conjugate to these flelds are: 

entropy S conjugate to temperature T 

volume V conjugate to pressure P 

magnetization M conjugate to mag. field h 

Of course V and P may be swapped round: either can be viewed as an external 



More common thermodynamic observables are the specific heats at constant pres- 
sure and volume, respectively Cp and Cy', the bulk compressibility, K; and the 
energy density, e. 

Equilibrium for given fixed external fields is described by the minimum of the 
relevant thermodynamic potential: 




^(k) ~ k-^-\ 



field. 



E for fixed 

F for fixed 

$ for fixed 

W for fixed 



S,T 
T,V 
T,P 
S,P 
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A phase transition occurs at those values of the external fields for which one 
or more observables are singular. This singularity may be a discontinuity or a 

divergence. The transition is classified by the nature of the typical singularity 
that occur. Different phases of a system are separated by phase transitions. 
Broadly speaking phase transitions fall into two classes: 

(1) 1st order 

(a) Singularities are discontinuities. 

(b) Latent heat may be non-zero. 

(c) Bodies in two or more different phases may be in equilibrium at the 
transition point. E.g., 

(i) the domain structure of a ferromagnet; 

(ii) liquid-solid mixture in a binary alloy: the liquid is richer in one 
component than is the solid; 

(d) the symmetries of the phases on either side of a transition are unre- 
lated. 

(2) 2nd and higher order: continuous transitions 

(a) Singularities are divergences. An observable itself may be continuous 
or smooth at the transition point but a sufficiently high derivative 
with respect to an external field is divergent. C.f., in a ferromagnet at 
T = T, 

(b) There are no discontinuities in quantities which remain finite through 
the transition and hence the latent heat is zero. 

(c) There can be no mixture of phases at the transition point. 

(d) The symmetry of one phase, usually the low-T one, is a subgroup of 
the symmetry of the other. 



An order parameter, ^, distinguishes different phases in each of which it takes 
distinctly different values. Loosely a useful parameter is a collective or long-range 
coordinate on which the singular variables at the phase transition depend. 
In a ferromagnet the spontaneous magnetization at zero field, M(T), is such an 
order parameter, i.e., 

M(T)=^fim M(h,T) 



then |M(T)| = for T > T^, and |M(T)| > for T < T^, 
Note: 3(T - T^)^ will not do. 
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^ is not necessarily a scalar, but in general it is a tensor and is a field of the 
effective field theory which describes the interactions of the system on macro- 
scopic scales (i.e., scales much greater than the lattice spacing). The idea of 
such effective field theories is common to many areas of physics and is a natural 
product of renormalisation group strategies. 

Examples 

(1) The Ising model in 3 dimensions. 

The theory is defined on the sites of a 3D cubic lattice and the variable on 
the site labelled by n is an e Z2. There is a nearest neighbour interaction 
and an interaction with an external magnetic field, h. Then the energy is 
written as 

H = J^an(Tn+^ - h^an, 

n,fi n 

where /x is the lattice vector from a site to its nearest neighbour in the 
positive direction, i.e., 

//e (1,0,0) (0,1,0) (0,0,1). 

The order parameter is the magnetization, 

^ n 

where V is the number of sites in the lattice. Note, that whilst the are 
discrete, M is a continuous variable, — 1 < M < 1 in the limit y ^ 00. 

(2) H2O 

Look at the two phases of liquid and vapour. The order parameter is the 
density, p, which is large for the liquid phase relative to its value for the 
vapour phase. 

The properties of both systems, their similarities and differences are best exhib- 
ited by showing the various phase diagrams. 
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PHASE DIAGRAMS 



Ising model HjO 



(b) 




Tc T 



h = 0± ^-^c(T) = 0± 

M = M*(T) P = P*c(T) 

h>0 and h<0 



II is the chemical potential, and /i = ij,c{T) is the line of first order transitions in 
the (/X, T) plot for H2O. It corresponds to the line h = in the {h, T) plot for the 
Ising model. 
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(i) Approach along (a) gives a 1st order transition whilst approach along (b) 
through the critical point gives a 2nd order transition. 

(ii) The order parameters are: 



magnetization M 



density p 



The conjugate fields are: 



magnetic field h 



chemical potential, ji or 
pressure, P 



(iii) The behaviour near T = T^{t= ((T - Tj/Tj 



(a) t^O-, h = {)± 



t^O-, II- iic{T) = 0± 



M(T) ~ \tf 



p{T) - - W+ 



Clear symmetry in curve 



No obvious symmetry but 
experimentally /3+ = /3_ 



(b) t ^ 0+, h = Q 



t ^ 0+, 11 = Hc{T) 



Susceptibility 



^ Ko{T) 



Ko{T) is for ideal gas 



Then 



x{t) ~ \t\ 



-7 



Note that for the Ising model with t 0-, /i 0+ we find x{T) ~ \t\^'^' 
with 7/ = 7. It should be remarked, howver, that 7/ is not defined for all 
models. 
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(c) t = 0, h^0+ 



t = 0, /X - ^ 0+ 



P- Pc^ {P- l^c)'^ 



/3, 7, S are examples of critical exponents, 
(iv) Coexisting phases 





(a) States between the curves I and 11 are physical but metastable. They 
do not violate thermodynamic inequalities. In the PV plot this is 
equivalent to 




which means that the compressibility is positive. This inequality is 
derived from entropy being a maximum in equilibrium. Howver, these 
states are unstable against changing to the mixed system, e.g., domains 
in the Ising model (or ferromagnet), and liquid-gas mixture for water. 
The continuous curves shown are the {h, M) and (P, V) curves for a 
pure phase. E.g., the Van-der-Waals equation of state: 



a 



P+-] {V-h)=cT 



(b) The Maxwell construction gives the true equilibrium curve taking into 
account the formation of the mixed system. The mixture is of the two 
phases A and B. The rule for finding the interpolation is illustrated in 
the case of H2O: 



Pa = Pb, Pa = Pb ^ Pa - Pb = v dP = 0, 

Ja 
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where v = V/N, and N is the number of particles. This is the equal 
areas rule of Maxwell. 

What happens is that boundaries between phases or domains form. 
In each domain the magnetisation is oriented differently and so the 
bulk average magnetisation can be any value in the range —\M\ to 
|M|, where |M| is the magnetization of a pure domain. The walls do 
increase the energy of the system by Ae and there is also an increase 
in the entropy, AS. since there are many ways of realising the mixed 
state. However, the resultant change in the free-energy, AF = Ae — 
TAS , depends on the surface area of the walls and is negligible in 
the limit of very large volume. Of course, the actual way in which 
domains, or bubbles, form and move is very important (e.g. in the 
early universe, cosmic string formation etc.) but needs more analysis 
than the embodied in the Maxwell construction. 

2.1 The general structure of phase diagrams 

A thermodynamic space, F, is some region in an s-dimensional real vector space 
spanned by field variables yi, ■ ■ ■ ,ys (e.g., P, V, T, /i, . . .). InY there will be points 
of two, three, etc. phase coexistence (c.f. A and B in H2O plot above), together 
with critical points, multicritical points, critical end points, etc.. Q is the totality 
of such points. The phase diagram is the pair {Y,Q). 

Points of a given type lie in a smooth manifold, M, say. The codimension, k, 
of these points is defined by 

K = dim(F) - dim(M). 

E.g., two-phase points have k = 1; critical points (points that terminate two 
phase lines) have k = 2. 

There do not exist any known rules for constructing geometrically all acceptable 
phase diagrams, {Y,Q): we cannot construct all the phase diagrams which could 
occur naturally. 

2.1.1 Structures in a phase diagram: a description of Q 

I assume that there are C components in the system, and hence there are (c+ 1) 
external fields: iJ,i, . . . , iJ,c,T. Then dim(y) = (c + 1). 

(a) Manifolds of multiphase coexistence. 

The Gibbs phase rule states that the coesistence of m phases in a system 
with C components has 

/ = c + 2 — m 

where / is the dimension of the manifold of m-phase coexistence. 

proof: dim(y) = (c + 1) and hence the manifold has codimension k = 
(c + 1 — /). But K = (m — 1) since k external fields must be tuned to bring 
about m-phase coexistence. 
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(b) An ordinary critical point has k = 2 and occurs when two coexisting 
phases become identical: 





two coexisting phases phases are identical: a = b 



critical point 

b 



(c) A critical end-point occurs with codimension k when two coexisting 
phases become identical in the presence of (k — 2) other phases: 





a = b coexisting with c 

K = 3 



a = b coexisting with c,d 

K = 4 



An example of all these structures in a three dimensional phase diagram is 
shown below. 
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first order surface = 
2-phase coexistence 




thus the critical end-point terminates a line of critical points and also ter- 
minates a line of triple points. 

Note that a phase diagram can often only be properly understood if plotted 
in the space of all relevant parameters. E.g., the Gibbs rule might seem to 
be violated since too many phases are coexisting at one point. However, if 
the space is enlarged in dimension this will be seen as a special case which 
only occurs for a particular cross-section of the enlarged space; 



a 



forbidden for c = 1 
but allowed for c = 2. 



(d) A tricritical point has k = 4. Its nature is most easily seen first in three 
dimensions. This is already a special case since we can only be sure it 
will appear in four dimensions. We suppose we have taken the appropriate 
cross-section of the 4D space, this often occurs naturally since some of the 
parameters are naturally set to the special values necessary to show up the 
tricritical point: e.g., by symmetry considerations. 




The hatched surfaces are Ist-order surfaces: surfaces of two phase coexis- 
tence. Thus the 1st order hne in 2D is really a line of triple points (three 
phase coexistence) in higher dimensions. 

A less special 2D cross-section of the same model will be: 



2 THE PHENOMENOLOGY OF PHASE TRANSITIONS 



13 




Here we have set h = (see 3D plot) but have changed the value of a fourth 
parameter u. 

(e) A bicritical point is a critical point at which two critical lines terminate. 
A typical phase diagram is: A model which has a phase diagram like this 



g 




bicritical point 



is given by the Hamiltonian 

<ij> ^ j ^ ^ ^ 

< ij > means nearest neighbour pairs, i.e., it labels the hnks on the lattice. 
Sj is a vector at the i-th site with |Sj| = 1. 

For low T thermal fluctuations can be ignored and it is safe to just flnd the 
configuration (i.e., set of values) of spins {Sj} which minimizes H. Since 
J > the first term causes the spins to align with each other to give 
ferromagnetic ordering. 



g<0 



Ordering is prefered along z-axis. This is phase B. 
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g > Ordering is prefered in the xy-plane ± z-axis. 

This is phase A. 
g = Neither A nor B prefered: two-phase coexistence. 



For high T ordering is absent: it is destroyed by thermal fluctuations. For 
g ^ as T increases we must have a second order transition from ferro- 
magnetic to paramagnetic phases. The surfaces to the low-T side of the 
Hues of critical points are flrst-order surfaces. This can be seen by impos- 
ing a magnetic held h on the system with components both || and _L to 
the z-axis and adding the magnitude, h, of h (including sign) as a third 
orthogonal axis to generate a 3D phase plot of which our 2D plot is the 
h = cross-section. Then as h changes sign the magnetisation, M, changes 
discontinuously at h = 0. This occurs as the surfaces in our 2D phase dia- 
gram on the low-T side of the critical lines are punctured, and hence they 
are in fact flrst-order surfaces. Of course, because the order parameter is a 
vector the possible patterns of behaviour and the competition between the 
effects of the terms governed by the coupling, g, and by h is, in general, 
complicated. An r-critical point is where r critical lines terminate. 

(f) A critical point of rj-th order has k = n + 2 and is complicated. 

3 Landau- Ginsberg theory and mean field the- 
ory 

The Landau-Ginsberg theory is a phenomenological theory describing all types 
of phase transition which can be derived from the more complete theory. It is 
a classical approach which breaks down in its simple form for low dimensions. 
However, it can be used for developing the structure of phase transitions and 
phase diagrams. Landau theory gives the correct prediction for critical indices in 
dimensions d > dc, where dc is a critical dimension which is different for different 
kinds of critical point. E.g., for an ordinary critical point dc = 4, and for a 
tricritical point dc = 3. 

Mean held theory is a method of analysing systems in which the site variable 
(spin etc.) is assumed to interact with the mean field of the neighbours with 
which it interacts. In a spin model each of the neighbouring spins has the value 
of the mean magnetisation per spin, m. The problem now reduces to that of 
a single spin in an external held and can easily be solved. By demanding that 
the mean value of the spin in question is m the solution yields a non-linear 
equation expressing this assumption of self-consistency and from which m can be 
calculated as a function of T. The approximation of the method is that it ignores 
fluctuations in the spins about their mean. It will turn out that Landau theory 
suffers from the same deflciency as we shall demonstrate. Mean held theory and 
Landau theory give the same, classical, predictions for critical exponents. 
We shall consider the following example. 
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Let the order parameter be M, and expand the free energy, A, as 
A = Ao + \a2M^ + \a^M^ + \a^M^ + . . . , 

^ rr U 

with ^ oc (T - Tc). 

There are no terms with odd powers of M. These can be present in principle 
but can be consistently excluded by symmetry considerations if the microscopic 
Hamiltonian is invariant under M — >• — M. If odd powers of M are present then 
generally the theory has only first order transitions, although higher order transi- 
tions cannot be totally excluded. Tc is a complicated function of the couplings in 
the original, microscopic, Hamiltonian as are the other coefficients, A2n- It is an 
assumption that A^ is analytic in T: an assumption that can only be plausibly 
justified under certain circumstances. This assumption as well as others is wrong 
if the dimension is low enough. 



Equilibrium is given by minimising A: 

dA 



dM 



0. 



The observable value of the order parameter, M{T), is the solution of this equa- 
tion. Then 



\M{T)\ = 



Thus as T ^ 



2/1 AeA2 



2 Ai 



\M{T)\ ~ |T-Te|2 ^/3 = 
We can rewrite the expression for M{T) as 



M(T) = 



A, 



where x = 



A.A 



6^2 



42 

^4 



If we assign to M a dimension of (—1) and dimension dto A then the coefficients 
A2n have dimension {d + 2n). Thus in this artifical dimensional analysis we find 
that X is dimensionless and that the critical exponent is predicted on dimensional 
grounds. The analysis above is only possible if A^ > 0, in which case A^ only 
occurs in the correction terms. If ^4 < then we require Aq > to stabilize 
the calculation and the results are different (see below). In the former case since 
only A2 and are important the critical exponent follows uniquely from the 
dimensional argument. 

A4 > 

If a field h is applied then the symmetry is broken and 

A = Ao-hM + -A2M^ + -A^M^ 
u '2 4 
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At T = Tc {A2 = 0) the condition for equilibrium is 

-h + AiM^ = Mr^h^ 5 = 3. 

For T > we have {t = ^^^) 

-h + a2TctM + A^M^ = 0. 
Then the susceptibility is given by 



X = 



dh y^^Q a2Tc 



1 T 

^ 7 = -i. 



The curve of ±M(T) vs T is a parabola: For T < we find 



. M 



parabola 




h= -a2|T-Te|M + ^M^ 

where in general A4 is a function of T which does not vanish at T = Tc. Thus 
the equation of state has the form 

h = -a{T)M + I3{T)M^ + . . . 

Then 




M M 



h>0 



h<0 
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If h is tuned from positive to negative the true minimum describing equilibrium 
changes from one to the other, but it takes time to re-establish th equilibrium 

especially if the intervening barrier is high. Consequently, the system can be in a 
metastable state corresponding to the local but not global minimum. This is the 
phenomenon of hysteresis. What happens is clear from the equation of state: 



(i) For hi there is unique minimum and the state is stable. 

(ii) For /i2 there are two minima (a) and (c) and a maximum (b). State (c) 
is stable and (a) is metastable, but (b) is unstable corresponding to a 
maximum of A and thermodynamic inequalities are violated here. 

(iii) Follow what happens as h decreases from hi through h2 and h^ and eventu- 
ally becomes large and negative: M varies smoothly as a function of h and 
metastability (supercooling) occurs, then at (d) the state becomes unstable 
and any fluctuation precipitates the change to the true equilibrium state. 

Plots of ^ vs M for the different h values make this interpretation clear: 



The Maxwell construction tells us that the stable state is not most generally 
characterised by a constant M. Locally M is a constant but it can change 
globally giving a domain structure. The Maxwell construction corresponds to 
the situation when h = and there are two degenerate minima associated with 
the two different domains that can co-exist. Thus at h = any value of M 
between these minima is possible and corresponds to an appropriate mixture of 
domains. The above analysis relies on the smoothness and differentiability of all 
functions A, h, M and can never directly address the mixed-phase system. 



Ah 





h, 







«0 
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A4 < 

As T decreases A{M) behaves qualitatively differently depending on whether 
^4 > or ^4 < 0: 




second-order transition at T = Tc first-order phase transition at T(, , T„> 



Hence the system passes from a second-order transition to a first-order transition 
as A4 changes sign and becomes negative. 
The stationary points are at M = and at 



-A^ ± [aI - 4^2^) 



/2Ae = Mi 



The + sign gives the minima and the — sign the maxima. 

To is determined by A{M) = having a double root at M = ±M+ (note that ^0 
is set to zero so that ^(0) = is the minimum for T > Tq). The solution is 



and at the transition 



V / + 4 4„ 



6 

Thus the point T = Tg, ^4 = separates the first-order line from the second-order 
line: this is a tricritical point. To see the tricritical point these two parameters 
have to take these special values and this requires tuning two external fields in 
the phase diagram. 

In the space of physical fields, denoted by T and g (e.g., g can be identified with 
a chemical potential controlling the relative abundances in a two component 
system), the phase diagram has the form: 
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TCP 



\ Tc(g) 



T 



Alternaltively the phase diagram in {A2, A4) space takes the form: vskip 5 truemm 



(i) Both trajectories (a) and (b) exhibit a second-order transition. 

(ii) Trajectory (a) passes through the TCP and lies entirely within the tricritical 
region. The transition is characterised entirely by the properties of the TCP, 
and all critical exponents are tricritical ones. 

(iii) Trajectory (b) exhibits an ordinary second order transition. However, it 
starts in the tricritical region and so initially the divergence of the relevant 
quantities is controlled by the TCP. Eventually it passes into the critical 
region and the transition is characterised by the line of ordinary critical 
points and the critical exponents that are given above. 

In other words we only see a transition controlled by the TCP when we apporoach 
along a trajectory lying in the tricritical region. For trajectories that pass from 
one region to another we see a change in the critical behaiviour. This change is 
characterised by crossover exponents. 



critical 
region 



tricritical 
region 



second-order 
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At the TCP = 0, and so we have 



\M\ = 

New exponents can be defined: 

(To - TJ 
AM 



/3t 



with A = ^■ 



(-^4)* with i> =1: 

1 

2' 



|^4p" with ^„ 



where AM is the discontinuity in M across the first-order fine (^4 < 0). 
then for small ^4 we can write 



IMI = 



re = 



m(a;), 



2|^|2^| 

To see tricritical behaviour along a trajectory we clearly need x small, i.e., 

<4|^|A. 

This defines a parabola in the (^4, plane separating the tricritical from critical 

regions. This is shown on the figure. In the space of physical parameters it 
translates into a similar shaped curve defining the two regions controlled by the 
TCP and ordinary critical points respectively. 

The general theory of continuous phase transitions can be encoded in terms of 
scaling functions and relies on dimensional analysis together with some assump- 
tions about the behaviour of the scaling functions for small argument. If naive 
or engineering dimensions are used this is generally a recoding of Landau the- 
ory but is often used to describe the behaviour of the relevant thermodynamic 
variables as a function of the actual external fields and hence parametrises the 
experimental observations. 

Add a magnetic field, /i, with contribution to the free energy of —hM. Then we 
can always write 



A = 



\A, 



Al 




The point is that the equilibrium free energy, A, can always be written in terms 
of dimensionless ratios in this way. As before assign dimension (—1) to M and 
dimension d to A, and then A^ has dimension (d+ri). The above expression is then 
a general way of writing the dependence of A at equilibrium on the coefficients 
A2n in terms of a scaling function, F . Note that since A^ is always taken as 
positive it causes no problem in the denominators. 
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We now compare with a standard form parametrising A near the TCP: 

9 h \ 



A=\T-Tm?-'' F 



\T-T,{g)\r\T-%{g)\^j' 

where ^ oc A4, and thus measures the distance from the TCP along the tangent 
to the critical line at the TCP. Note that g has been substituted for the field g 
as the second independent external field: the critical line is thus parametrised as 
Tc{g). The TCP then is at position (0,Tt) in the T) plane where Tt = Tc(0). 
Labelling the critical exponents at the TCP by suffix, t, we clearly have 

1 ^ 1 A 5 

The following examples clarify the interpretation. 

(i) /i = 0, ^ = 0, T ^ Tt such that = a; is fixed. Then 

A = \T -Tt\^F{x,0) with F(0,0) finite. 

We see tricritical behaviour and since g ^ \T — Tt\'^'^ the trajectory fies in 
the tricritical region, (pt is the cros-over exponent. 

(ii) h^O, g fixed, T T^. 



1 3 



A= |T-Te|2 F 2 r,0 . 

V|T-T,|r J 

the argument of F is not under control and so we rearrange the expression: 

A = J ^ G ^,0 

9 \ 3 

where G{z,0) = zF{^,0) and G(0,0) is finite and non-zero. This property 
of G is an assumption in the general theory and could be violated. It does 
follow from the standard Landau analysis and hence if it turned out to be 
false in an experiment it would signal a breakdown of the Landau theory. 
The goal then would be to rescue the dimensional analysis approach by 
assigning values to the dimensions of the parameters different from the 
naive ones but which render the scaling functions F and G well behaved 
for small argument. 

In this case we find that A shows the normal critical behaviour associated 
with an ordinary critical point, namely 

An.\T- T/"° with a = 0. 
(iii) T = Tt, g^O, hj^O 



A = lim |T-T.p-"' F 



9 h 
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but A remains non-zero iov h ^ and we must be able to rewrite A in the 
limit as 

A = h~ J 



i.e., F must have a singularity that cancels the external factor of |T— T^p 
Thus 

A = hlj(A 

If we assume J(0) is non-zero and finite then at the TCP 

A and M = — => St = 5. 

oh 

This result can be derived much more directly from the Landau theory: at 
T = Tt both A2 and A^ are zero and so we can write 



A = -hM + AqM^. 



Thus we find that 



dA 

m 







However, it is important to see how the parametrisation in terms of scaling 
functions works, and, as has already been remarked, this form of parametri- 
sation is more general than the naive Landau theory: the assumption that 
F(0,0) etc. are non-zero and/or finite breaks down when Landau theory 
ceases to be valid. 

A summary of the critical indices is 





An. |t|2-° 






M 


cross-over 




a 




7 


5 




CP 





1 
2 


-1 


3 


1 
2 


TCP 


1 
2 


1 

4 


-1 


5 


1 
2 



(To-Te) ~ \A,\^ ^ = ^ 



AM 



To finish this section we look in general at the Landau model with one parameter. 
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Let the thermodynamic potential be ^ and the order parameter x. Then 

^ = OO + CLl^ + CL2X^ + ■ ■ ■ x^''. 

We assume that the coefficient of the highest order term is 1. This means that 
higher order terms will not affect the phase structure since the coefficient of this 
term will not be allowed to vanish at any point in that part of the phase space 
in which we are interested. In addition this highest power is even. This ensures 
stability, i.e., has a global minimum characterising equilibrium. 
The Oj are thermodynamic field variables which span the phase space, Y, and 
for a given choice of the the stable thermodynamic state is identified with the 
value of X which minimises ^. Choose gq such that this minimum is always at 
^ = 0: this simply fixes the origin of ^. Then we can write 



where Vj the bj are real and the dj are real and non-negative. However, at least 
one of the dj is zero since ao has been chosen appropriately. Then ^ can be 
represented by the configuration of a set of points {bj, dj) in the upper half of the 
(6, d) plane. Consider the following examples. 

(i) m-phase coexistence 

The minimum of ^ occurs for m distinct values of x. These are identified 
with the m coexisting phases and it corresponds to m of the dj being zero. 
In general this requires (m — 1) constraints on the dj and hence the man- 
ifold of m-phase coexistence has codimension k = m — 1. For m = 2 the 
configuration takes the form: 



^ b 



(ii) A critical point 

A critical point occurs if two of the dj are zero and the corresponding bj 
coincide: 
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> b 



(•) = two coinciding points 



For this configuration 



* = (x-bY + a5{x-bf + .... 



There is no quadratic term (rr — 6)^ and hence have a second order transition. 
This is a more general situation than before where we imposed symmetry 
under x — >■ —x for simphcity. 

The additional constraint on the bj means that k is one more than for two 
phase coexistence, i.e., k = 2. 

Remarks 

(a) All points of m-phase coexistence can be continuously transformed into 
each other: inspect diagrams of the kind shown in (i) above. Hence 
they all lie on one manifold. 

(b) There is only one manifold of critical points. This is clear from (ii) 



above. 



(iii) Critical end-points 



4 d 




b 



or alternatively 
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4 d 



* b 



It is not possible to pass continuously between these two diagrams and 
hence there are two distinct manifolds of critical end-points separated by a 

(iv) Tricritical point 



This point has k = 4: two constraints occur because there are three coexist- 
ing phases (i.e., m = 3 =^ (m — 1) constraints), and two further constraints 
arise because of the requirement that the relevant bj coincide. Howver, k is 
just the number of such constraints. 

(v) Multicritical points of order r 

r of the dj are zero and the corresponding bj coincide, thus 



clearly we can determine how the manifolds fit together by following the 
fiow of points in (6, d) space as the vary. 



4 d 




K = 



(r - 1) + (r - 1) = 2r - 2. 
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4 The Partition Function and Field Theory 

The partition function Z is defined by 

where ^(x) is a generic field degree of freedom and /3 = 1/kT. H{<f)) is the 
Hamiltonian given by 



H{(t)) = / dx 7^((/)(x)) 



where ?^(0(x)) is the Hamiltonian density.A is the large momentum cut-off which, 
for a lattice of spacing a, is A = 2w/a. In this case the integral will be a sum over 
all sites of a discrete Hamiltonian density. The crucial point is that there will, 
in general, be a cut-off of some kind. 

We shall assume that the coefficients in H depend only on the volume, V, of 
the system. Then, for a given temperature T and volume V of a subsystem the 
equilibrium probability density for finding the subsystem with field configuration 
^(x) is 

Then the entropy S is given by 

s = -k j {d<i)}pmog^m 

= -k j {(i(/'}^e-^^('^^ {-I3H -log Z) 



= k{l3U + log Z ) 



where 



Thus 



and hence 



^ ~ j {d^}H{(j))e = internal energy. 



A;Tlog Z = -U + TS = -VF, 



F is the thermodynamic potential appropriate for T and V as independent vari- 
ables. 



We are interested in the macroscopic properties of the system and so we re-express 
Z in terms of a macroscopic variable (/)(x). This is our "guess" for the order 
parameter, but it could be that it will not reveal all the possible phases of the 
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system. We are trying to pick out the collective coordinates on which the long- 
range/low- momentum physics depends. For a scalar field system this is not too 
difficult but when the field has internal degrees of freedom it is much harder. For 
example, both the phenomena of superconductivity and He3 superfluidity have 
spin-i fundamental fields which pair to form spin-0 bosons which condense. BUT 
in the superconductor this pairing is in the state whereas in He3 it is in the 
^Pq state. Hence different choices must be made in these two cases. 
We can choose 

where v is a volume centred at x with v = ^ a^. 
Alternatively 

y|p|<A<A [Zir)^ 

Here 0(x) is composed of the low- momentum degrees of freedom only : A = 27r/L. 
The cut-off is now A/L and hence it is possible to define the field theory for (j) on 
a lattice of spacing L. Then 

Suppose, for example, 

7^((^(x)) = i(V<^)2 + ^m^<^2 + i^o*/*' . . . + lj9nr. 

where n is taken as even and where itiq, gn are the bare or microscopic coupling 
constants; i.e., they parameterize field interactions on a lattice of spacing a = 
27r/A. H{L, (j)) is much more appropriate to our needs since the unit of length is 
now L ^ a. Write 

nm 4>) = T)iy4>f + ^-m\L, T)4? + i^(L, T)4>' + . . . + (V^)^ + . . . 

Note that the kinetic term no longer has unit coefficient. Generally at this stage 
it is convenient to rescale cj) to make this particular coefficient unity: 

Z5(L,T)<^. 

However, we will not do this rescaling for the moment. 
Typically as L becomes large enough 
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Hence the effective Hamiltonian is insensitive to and so ^ can be treated as 
a constant. Then 



H = lm2(cx),T)^2 + ^^(cx),T)(^^ + ... 

We have integrated out over all scales and produced renormalized coupling 

constants and an effective Hamiltonian depending on 4>- 

Now 

1 



^ g-VS(xo) 



e-^s(xo) 



27r \ ^ ' 

y5" 



1 - 



where ^'(xo) = 
Hence 

'logy^ 



F = H{oc,(t)o) + O 



V 



where 0o is the global minimum of H{oo,(j)). This is Landau's method with 
H{oo, (j)) as the free energy, function BUT this procedure requires the limit 
L — )■ oo to be under control. In particular Landau assumes that m^(oo,T) is 
analytic in T. Then all that the integration over scales has done is to fix the 
value of Tg. This assumption is wrong for D < Dc and Landau's method fails. 
Note that (j)Q is the analogue of the magnetization: it is the order parameter 
and can be measured in an experiment. 

Since it is the limit L — t- oo which causes the trouble let us investigate calculating 
F with L finite and identify the source of the problem. 
Let S{4>) = l3H{L, (j)) and define by 



= 0. 



Then 



where 



^ fdydz (-^f^rA ^(y)^(z). 
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Assuming 4>o is a constant independent of position and taking logs we find 

F = n{Ljo) + log l^e-^^"^' + - . 

In H the mass and couplings are analytic in T for fixed L. (This is certainly 
plausible and is a reasonable assumption). The integral is calculated using the 
loop expansion. Landau assumes that these loops just change Tc and that for 
L large enough they may be neglected keeping only the first term in F above. 
To illustrate how problems occur consider the effect of an external constant field 
J. 

The induced magnetization is 
and the susceptibility is 

That is, the integral over the connected two-point function 

< (/.(0)(/.(x) >e=< (/.(0)(/.(x) > - < (/.(O) >< (/.(x) > . 

The second term arises from differentiation of Zj in the denominator. 
We now define the smoothed field (j){L, x) including a renormalization to keep the 
coefficient of |(V^)^ in H{L,4>) to be unity. Of course, we must remember this 
rescaling if we re-express results in terms of the actual magnetization. Then 

_i 

X = ^Z,{L) J < 0(O)0(x) >e . 

We define 

Tl is the truncated two-point function. 

[ In general, the n-point function G„(pi, . . . ,p„) is defined by 

G„(pi,...,pJ = JdXi...dXri < (/'(Xi)(/l(x2) . . . (/>(xj >c e^iPi-Xi ^ 

and the truncated n-point function is defined by 

Gn(Pl, • • • ,Pri) 



r„(pi,...,Pr*) 



G(pi)G(p2)...G(p„) 
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Gn contains the denominator of the RHS as a factor and r„ is often a more useful 
quantity with which to work.] 
Then we have 

From now on we will absorb a factor of into ^, i.e., replace 

(f) -)■ ■ 
then ^ , g ^ (3g etc. and 

^ ML) ■ 

To calculate f ^ we need to use perturbation theory. 

4.1 The Perturbation Expansion 

We start with the definition 



< 

Let 



where can be negative. The situation < corresponds to spontaneous 
breaking of the symmetry —(j) in the bare Hamiltonian. However, we shall 
see that the effective mass is positive once loop corrections are included. 
In order to be able to perform the loop expansion we write 



and choose /x^ > 0. 

Then with 



we have 



^Jdct) m^i^) [l - / d^lim' - + i^c/)^ + . . . 
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This consists of a sequence of gaussian integrals which are best done in momentum 
space, they take the form 

J d4> e-J'^y^^ mQ{y-^)m (/,(0)(/>(x) 
= ^(detQ)"^ Q-\x) 

and 

jd(t)e-I "^^"^ (/'(0)(/>(x) j dz g<i)\7) 
= gB{detQy^ J dzQ-\z)Q-\y. - z)Q-\0) . 

The denominator cancels the (det Q)^ leaving a calculable combinatorial coeffi- 
cient. In momentum space we have the diagram expansion as follows 

G(p) = 1 c?x < (/)(0)(/)(x) e-'P^ 
or 

r(p) = Q-\v) + (m^ + <3-^(q) + . . . 

0^(p) is the bare Feynman propagator: 0^(p) = (P'^ + A*^)^- 
Then 



G(p) 



K 

P P 



r(p) = ( — • — )-• = ( + 



Cg 



+ Cg 




p p 

q 




The expansion shown for V is exact to 0{g) but corresponds to a sum over selected 
1-particle irreducible graphs in G. this can be seen as follows. 



r g + 5m2 + E ' 
where E represents 1-particle irreducible (IPI) diagrams. This means that 
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they cannot be separated into two pieces by cutting one line only. E.g., Expand 




n 



IPI IPR 



the above expression to give 



Diagrammatically we have 



To 0{g) : 



+ — X X — + — H @)— 

+ -@ H + -@ @- 



( • )' = ( )■' + + -@- 

T 





IPI 



Note that F is truncated because it contains no propagator factors, Q ^, for 
external legs. 

The expansion for r(p) can be improved by summing over all bubble insertions 
on the internal loop. Then 



r(p)= ( — • — r = P+m^ + Cg 







+ 
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The Feynman rules for constructing the diagrams are 



1 



for each propagator line 



S (27r)^ ^^^^ closed loop 

S^{YlPi) momentum conservation at each vertex 

factor C = ^ 

Where S is the number of symmetry operations that leave the graph unchanged, 
for example 



S = 2 f ) S = 6 





Note that disconnected contributions to the numerator, e.g.. 





are cancelled by terms in the denominator. Hence, only get connected graphs 
for expansion of G. For our smoothed theory we have 

= m'^{L, T). 

Remember that is not necessarily positive. Then 

^^^Aj^^T)_ 
p2 + m2(oo,T) ■ 

Here m^(oo,T) = ^^^(T) and is positive. The loop correction renormalizes 
additively and the resulting effective mass is positive. Now 

, fi(0) m2(cx),T) 

X = - 



Zi{L) Zi{L)Zl{0,T) 
X^^ vanishes at T = Tq which is given by the vanishing of the renormalized mass: 

lim m^(oo, T) = 0. 



Note that in general lim^^r^ Zl{0,T) = as well but that the numerator domi- 
nates. In fact for our present discussion 

Zi = l + 0(/) , 
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and hence to 0{g) Zl = 1. 
Then 



m (oo, T) = m (L, T) + -g —— — — — — + . . . . 

2 J (ztt)^ + m^(oo, T) 

Note that the bubble summation removed all reference to /i^. In fact we can 
see that by choosing /i^ = m^(oo,T) in the first place we automatically get the 
improved formula from the one-loop expression, the coefficient of | comes from 
the symmetry factor as described above. Using m^(oo,Tc) = we write 



2 d p 1 1 

m^(oo,T) = m^{L,T)-m^{L,Tc) + -g J j—r^ 



(27r)^ [p2 + m2(oo,T) p2 

= At--gmioo,T)J ^ (p2 + ^^(^^ r))p2 +••• 

Where t = {T -T^)/T^. 
We consider two cases: 

D > A The integral is dominated by large p and is finite in the limit 0. 
Then 

T) = At + Bm'^{oo,T). 
This is consistent with the Landau behaviour 

m2(oo,r) = Ct 

D < A The integral is sensitive to low p and is Infra-Red divergent unless > 
0: 

m2(oo, T) = At + Bm^-'\oo, T) 

This is clearly inconsistent with the Landau assumption: the integral is 
important. IR divergences invalidate Landau's assumption. 

D = A K marginal case with IR logarithmic corrections 

m^(oo,T) = At + Brm? logm^. 
Get Landau behaviour modified by logs. 

Hence Dc = A and for D < 4 IR effects destroy mean-field predictions. 
This theory applies to an ordinary critical point in Landau's theory: the coeffi- 
cient of (jp vanishes. If, for some reason, the <j)^ were absent and the interaction 
started at the g(j)^ term then the loop contribution would be 



/•¥ d^p 1 /■¥ d^q 1 

J (27r)^ p2 + m2(oo,T) J (2^ q2 + 



(oo,T) 
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Subtracting off tlie value at = and isolating the most IR divergent integral 
gives 

J (27r)^y (27r)^ p2(p2 + m2(cx), T))(q2 + m2(oo, T)) ' 

There are no IR problems for 2D — 6 > 0, i.e., Z) > 3. Hence in this case 
Dc = 3. In Landau's theory this situation corresponds to a tricritical point: the 
coefficients of and vanish. 
In general, for interaction we find 

{n - l)Dc - 2n = 0, 

or 



This interaction is relevant only for critical points of order n. For D = 2 all 
critical points show anomalous behaviour but for D = 3 only critical and tricritical 
points are anomalous, and even then the corrections to mean-field predictions at 
the tricritical point are only logarithmic. 

It is important to understand heuristically or phenomenologically what is hap- 
pening before actually doing detailed calculations. First we define the correlation 
length, ^, a bit more explicitly than before. 
For X ^ a 

< (/'(0)(/.(x) e-l"l/« . 
In fact if X = {xq, . . . ,xd-i) then 

J dxi... dxD-i < (/'(0)(/>(x) Ce"^°/« , 

as xq becomes large. This is how we might calculate ^ on a computer. 
[ Note 

J p^ + 

] 

Now, <j){L) averages over blocks of size L^, but if D < M{L) = Zl {L)^{L) 
cannot be identified with the magnetization of the bulk volume, V, since there 
are still fiuctuations in the system of wavelengths greater than L which cause 
these blocks to interact. In other words a block of size is not big enough to 
be a good model of a larger system. Only for L ^ will this be true. For D < Dc 
there are important fiuctuations on all scales L < ^, and only for D > Dc are 
these fiuctuations suppressed and Landau's theory works. 
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(1) In the theory with L = ^ we do not expect significant loop corrections. 
Hence 

< <^(0)<^(x) >e~ e-™(«'^)l''l , 
which identifies = m(^,T). 

If we postulate that m(^, T) ~ • (T - Tj, then we find 

^ - ^ = ■ 

l + a 

Also 

;^-i = r(p = 0)-^m2(e,T) . 
Then if Zi{^) ~ = tr-p we get 

(2) In general we have 

G(x) = < (/)(0)(/)(x) >e 

= ^(" \|x|eF-i)P a«^«|x|, 

= ^'(«)^^' = (i:> - 1 - p)/2 , a<|x|<^, 

where A is the lattice spacing, i.e., A = 27r/a is the UV cut- off. In all we 
do A is finite but A ^ m and for all momenta p of interest A |p| (i.e., 
a «C |x|). Hence we are always working in the limit of very large cut-off. 

(a) It is important to realize that ^ sets the scale and not a, otherwise the 
correlators would be badly behaved as ^/a oo and/or |x|/a — )■ oo. 
This is the basic assumption of the scaling hypothesis. 

(b) "Naive" dimensional analysis must work. 



= 0^[(/)] = ip-2) . 



Where [] signifies the dimension of the enclosed quantity in units of 
momentum. Hence 

[<(p(p>] = D-2 . 
Thus in the above expression for G(x) we must have 

C(a) ~ a'^ , r/ = - ^{D - 2) . 

r] is the anomalous dimension of the field <p, and A<^ is the scaling 
dimension of <f). 



4 THE PARTITION FUNCTION AND FIELD THEORY 



37 



Now 



X = j ^(^) J dx 



e 



X 



Comparing with the alternative expression for the behaviour of ^ in (1) 
above we find the scahng relation 

a = u{D - 2A^) . 

(3) The results of (2) can be seen in a different way. We have the general 
parameterization 

G{P)= ^'"'^ 



We were concerned with p = earlier (x = G(p = 0)), but now instead 
consider A |p| m. Because Z is dimensionless we have 

Z r-. (ap)-"' f (^j^ p:$>m = r\ 
where p = p. Also from above in (1) we have 



Z ~ (am) ''gi—] p<^m = ^ 



This is because Z ~ Zi{^) ~ for p = 0. 

Note that / and g do not depend explicitly on {ap) or (am) for the reasons 
outlined in (2a) above: because of the scaling hypothesis. Then 

(a) Z{p = 0, m) is non-zero and finite for m ^ 0. 
Z{p, m = 0) is non-zero and finite for p ^ 0. 

(b) the scaling hypothesis asserts that m (i.e. ^) sets the scale and not A 
(i.e. not a). Hence the above scaling forms must be vahd with /(O) 
and g{{)) non-zero and finite. It then follows that 



By comparing powers of a we also get 

P = -277. 

Substituting p in terms of a from part (1) gives the scaling relation derived 
in (2). 



4 THE PARTITION FUNCTION AND FIELD THEORY 



38 



The important thing to notice in all of this is that each quantity in which we are 
interested depends on a in a simple and special way. Namely, a appears raised 
to a power as a multiplicative factor only: a does not appear as part of an 
argument of any of the functions. These function depend only on the long range 
observables. Since a has dimensions it follows that whilst "naive" dimensional 
analysis holds when including the factor depending on a, the dimension of the 
function it multiplies is not constrained. This means that if we concentrate on the 
dimensions of this function (e.g., f or g above) the contributions to its dimension 
from the observable quantities can add up to something other than the "naive" 
dimension, i.e. the anomalous dimension. Since, in any given system a is 
fixed the anomalous dimension is the most natural dimension with which to be 
concerned. After all, we can vary ^ by varying T but a is, of course, unchanged. It 
is the object of most analyses of critical phenomena to calculate these anomalous 
dimensions. 



4.2 The Effective Potential and the Legendre Transform 

Consider the partition function in the presence of an external source J: 

Z{J) = jd(j) e-^^f^^-'l' . 

Denote J) = - J(/>. 
Then 



where 



and 



^ ) {(j)^ J) = when cj) = (j){J) = (/)(x; J), 



^(y) = Hy)-Hy;J)- 

Space arguments will be suppressed wherever possible. 
Thus 

\ogZ{J) = -5((/)(J), J) + log J d^j e-^^i'''^' 

The integral generates loop corrections: this is the loop expansion. It is essen- 
tially the same as the perturbation expansion we saw before except that here the 
expansion is about the minimum of S'(^, J). Then we have terms of the form 

/ e-'^''^'^' e-hs'^'^-')^' - = jdi; e-'^''^^' (l - ^S^'\j)^' - . . .) 

= A[dei{Sf)Y\l + ...). 
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We get a set of gaussian integrals which generate the loops. For example, the 
1-loop contribution to log Z{J) is 

-i [det(sf )] . 

In general 

^ ^/ ^ f S^^og Z{J) \ 



j=o 



(See problem sheet 2a for an explicit example of the 1-loop calculation.) 

The tree diagrams are generated by —S{<f>{J), J). Once we understand how this 

works we can verify the preceding statements. Consider the approximation 

W{J) = logZ{J) ^ Wt{J) = -S{(p{J), J) . 

4>{J) satisfies 



= . 

</,=</,( J) 



This is the field equation and 4>{J) is the classical solution. Then 

SWt{J) _ _JS_5m ^ MJ) ^ ^fn 
SJ ~ 5(i){J) SJ ^ SJ ^ ^^'^> 

= m- 

This field equation is explicitly of the form 

(- V' + m2)(/.(x; J) + ^^(/.^(x; J) = J(x) . 

Clearly all derivatives of W with respect to J can be generated from the previous 
equation once this classical equation has been solved. In perturbation theory 
we have, in momentum space 



p2 + m? 6(p^ + m^) 

Or, diagrammatically 



e 



jpx 



■^J - I □ (l)' 



Iterating this equation gives 
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+ — 

12 




All momenta are integrated with momentum conserved at each vertex. This is just 
the tree expansion of W{J). It is generated by the classical field solution, 
(j){J): i.e., there are no h factors. 
E.g., 



( '-^^^ ^ =o5^ 

V5J(p)5J(qi)5J(q2)5J(q3) 



= gS (p + Qi + q2 + Qa) 

j=o 



This is the tree contribution to the 4-point vertex. 
Now, two important points 

(1) If Wt{J) is the tree approximation to W{J) as outhned above, then Wt{J) 
is the Legendre transform of S{<f>). That is, it follows directly from the 
above discussion that 

Wt{J) = max{J(l) - S{(t))) , 
<t> 

which is the Legendre transform of S<f>). 
Then define 

SH{(t)) = max(J(/) - Wt{J)) . 
Clearly Sh{4>) is not a function of J, i.e.. 



dJ 



= . 



If >S'(^) is convex then Sh{4>) = S{<p), otherwise Sh{4>) is the convex hull 
of 
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This subtlety has a lot to do with the Maxwell construction and domains 
associated with a first order transition. It applies to a complete study of 
the symmetry broken phase but we shall not pursue it further here. 

(2) For (j) = M, a constant, S{M) was the effective potential which is min- 
imized to give Landau's approximation to the free energy. Remember the 
form of S{M) with an assumption of analyticity of the coefficients in T is 
what was needed: loops just modify Tg. Hence Landau's method is clas- 
sical: it corresponds in form to tree diagrams. In general, is the 
effective action of the classical or Landau theory. Note that <S(^) is ex- 
actly the action which appears in the expression for Z 

That is, this S{(f)) is the Legendre transform of Wt{J)- There are some 
niceties concerning whether one should use S or Sh here, but we will not 
discuss this further. 



W{J) generates the connected diagrams. Then schematically it must have the 
diagrammatic form 




The vertices of the trees are the IPI diagrams of the theory, i.e., they contain 
sums over all loop contributions but that cannot be separated into two disjoint 

pieces by cutting one line only. 

We can generate these tree diagrams with an appropriately chosen effective 
action. This is the effective action of the full theory. It is the Legendre 
transform oiW{J): 

Se{(I)) = max{J(f)-W{J)) , 
n n. J 

Hence SE{(f>) generates the IPI diagrams of the theory. 

Se{4>) is now treated as a classical action since it is all that is needed to generate 
the whole theory by its tree expansion. The equilibrium free energy is 

p _ -w{o) 

13V ' 

and 

W{{)) = max{J(t>-SEmj=o, 

= max(-5£;(^)) . 

<t> 

Hence the free energy function F{<f)) is identified as 
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and equilibrium corresponds to minimizing F{(p) with respect to 4>- This is like 
the Landau method but we cannot assume that the effective coupling constants 
in Se are analytic in T: they are only for D > Dc- 
Note 

(1) Se{4>) is always convex and, in fact, is the correct effective action describ- 
ing the system for all phases. 

(2) If J is a constant function then the magnetization is given by 

M = 



dJ 

(/>(J) = M . 

In general, J(x) and M(x) = are conjugate variables. This is also 
shown by how they appear in the Legendre transformation. This is, of 
course, the same relationship that they have in thermodynamics. 



4.3 The Ising model 

The Hamiltonian is given by 

— = SiVijSj + JiSi . 

Si is that variable at the i-th site with 5, = ±1. Jj is the external field and Vij 
contains the couplings between neighbours: it need not be nearest neighbour. 
Introduce Xi E TZ and note the identity 



= Const • e^'^'^^^ . 



Then 



Si 

Now sum on Si 



s=±i 

Hence, ignoring all irrelevant multiplicative constants 



We calculate the free energy by steepest descent (i.e. by expanding about the 
minimum of the exponent as we have done before) and keeping only the first term 
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is equivalent to Landau theory: that is, to the tree approximation. Minimizing 
the action gives 

This defines Xi which minimizes the action. Thus in this approximation 
Wt{J) = log Z{J) 

= j,)y..i(x^.- j^.) + ^A{X,). 

i 

But 

Slog Z{J) 



SJ^ 

1 

2 



and using the defining equation for X above we get 

SA - 

M, = — (X,) = tanh(X,) . 
oXi 

(Remember, X and M are both functions of J.) 

To find the free energy function F{M) we take the Legendre transform of Wt{J) 

-I3NF{M) = max (J • M - Wt{J)) , 
where N is the number of sites. Then 

(3NF{M) = -VijMiMj + i^[(l + M,)log(l + MO + (1 - M,) log(l - M,)] . 

i 

Note that (3NF{M) in this approximation is the correct quantum action to be 
used in the full theory which includes all loops. This is because it defines the 
tree diagrams of the theory and hence gives all the bare vertices, the loops are 
built up in the usual way from these vertices. 
For constant magnetization Mj = M and M small 

F is the free energy per site and V = J2j Vqj. Then 

(3F{M) = (1 - V)M' + ^M' + ... 

Since V = V{j3) we have a second order transition at j3 = /3c given by V{j3c) = |. 
In the simplest case 

Vij = ^Kij^ 13c = , K = Koj- 

3 

Notes 
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(1) This is the mean-field approximation and we know that the critical expo- 
nents are analytic only correct for D > Dc- 

(2) The external field method together with the Legendre transform enabled us 
to calculate the free energy as a function of the relevant order parameter, 
M. Near the transition the Ising model can be described by a scalar field 
theory with the full theory based on the action F{(f)). F is the function de- 
fined above by F(M) for general (i.e. not constant) M. ( (f)is identified with 
M to make the correspondence clear). Thus we have found the action for a 
scalar field theory which is, in every way, equivalent to the Ising model near 
a continuous phase transition. This is the phenomenon of universality. 

Universality: many models which have disparate descriptions on the mi- 
croscopic scale exhibit the same critical properties with the same critical 
exponents. Near critical points these models are described by the same field 
theory and belong to the same universality class. Here we have found 
that the Ising model and scalar field theory belong to the same universality 
class. 

(3) Even for D > Dc loop contributions will give corrections to j3c. The 1-loop 
correction to the tree-level effective action is given by |detF^^)(^). This 
contribution is the subject of example sheet 2a. 



4.4 Calculation of the Critical Index 

This section is concerned with the calculation of the critical index i/ in the Ising 
model. We work with field theory in D dimensions. 

5(<^) = |dxi(V(/))2 + imV + ^#' +■■■ 
= 2 [g] = 4-D = e. 



We may imagine that we have obtained this action by integrating out the short 
wavelength modes. In general all types of interaction are present. Clearly ^ has 
been renormalized to give the canonical form to the kinetic term. Alternatively, 
we can accept that 5(0) has been generated from the Ising model by the trans- 
formation described in the previous section. We note for interest that the Ising 
model is identical to field theory in the limit vr? — >■ — oo, g ^ oo such that 
Glm^l/g = 1. 

Suppose that we have integrated over all momemta p > A and that the effective 
parameters are m^(A, T) and g{A,T). We shall also allow for multiplicative field 
renormalization constant Z{A,T). Now integrate over the next momentum slice. 
Write 

(/.(x) = (/.>(x) + (/.<(x) . 
0> contains contributions from momenta in the range A — 6A < p < A. 
0< contains contributions from momenta in the range < p < A — SA. 
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Then 



Z = J d(t>^d4>y e-^(^>+'^<^ . 
We do the integrals over ^> to get an action as a function of ^< only. 

S{ct>^+ci)y) = S{(j)^)+S{(j)y) + j d^[Vcj) 



+ 



12 



y(ix[2(/>^(/.>+3(/)^(/)^ + 2(/><(/>|] . 



The term quadratic in the fields vanishes identically. Since in p-space ^> and ^< 
have disjoint support 



j (ix(/)>(x)(/)<(x) = j c?p^<(p)^>(p) = 



Then 



As before we expand all exponentials of non-quadratic terms 



Z = 



12 



2 



The integrals over odd functions of ^< vanish and so will be omitted from now 
on. 

After the dcj)^ integrals have been done we are left with a polynomial in ^< 
and its derivatives. We gather these terms up and absorb them into 5(^<) by a 
redefinition or renormalization of the coupling constants. 

The important one-loop terms and their diagrammatic representation are 



-ijd^4>l4>l 




^^^ld^<\>l<\>\\dy<\>l^\ 




^'<(y) 



12 




<^^<(y) 



4 THE PARTITION FUNCTION AND FIELD THEORY 



47 



1-| is the propagator of ^> and has support in p-space only for 

A — SA < p < A. The last of the terms above is zero since ^< has disjoint p-space 
support from that of the intermediate propagator. 

We can now read off the renormalizations of m^(A,T) and g{A,T) to 1-loop 



m (A - SA, T) = m (A, T) + — - — / — -— 

2 27r 



/a-ja (27r)^ p2 + m2(A,T) ' 
In momentum space we have for the 4-point function 

^2 f d^'q 



where 



ASA {271)^ p2 + m2(A,T) (p + q)2 + m2(A,T) ' 
with the restriction A — 5A < |p + q| < A. 

We can expand r(q) = r(0) + ^^^(O) + . . . and then the operator generated is 



16 



r(o) J dx(t)\ + 4T'{o) j d^{4><y4><) 



Note that new operators are generated as well as ones we are interested in. 
It should be remarked at this point that the analysis of the higher dimension 
operators generated in this way is not easy and is better approached by other 
renormalization group techniques which will be touched upon later. However, 
there does exist work that has pushed the above style of analysis, i.e., momentum 
thinning, successfully to the study of the role of such operators. To 1-loop we do 
not encounter these difficulties. Also note that because the graph 




is independent of q we do not generate the operator (V^<)2 to 1-loop. More on 
this later. 

Then we find the equation for g 



g{A-SA,T) = g{A,T)-^g\A,T) 



d^p 1 



A-5A (27r)^ (p2 + m2(A,T))2 



We note that from the equation for the renormalization is positive. Hence 
if we want 

lim m2(A,re) = , 
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we require m^(A, T) < 0: the "bare" mass is negative. The expansion is still fine 
since it is only necessary that the gaussian exponent gives convergent integrals. 
This demands only that 

A^>-m^{A,T) VA . 

It can be verified in what follows that this inequality is always respected. 

We now rewrite the two renormalization equations in difi^erential form in terms 

of dimensionless quantities. 

u\b,T) = A-W(A,T) 
X{b,T) = A-'g{A,T), 

where b = log(Ao/A) and 6 = 4 — 0. Aq is the cut-off for the original theory. 
Then we have the evolution equations 

= 2u^ + 



db 2(27r)^ 1 + u 

dX 3^D 
= eA 



2 



db 2(27r)^ (1 + ^x2)2 ' 

is the surface area of a unit sphere in D dimensions: = 27r^. 
In order to integrate these equations the propagator for ^< must take the canon- 
ical form l/(p^ + m^). In general, the renormafized action, 5ij(^<) will be of the 
form 



- i^Oog Z)(V(/)<)^ + ^V(</'<V(/)<)^ + . . . 

Thus we must rescale the field to give the canonical quadratic part. Thus the 
new field for the next iteration of the procedure is 

0(x) = (1 - i^«)0<(x) . 

This defines the renormalized field <f) for finite b 

(/)(x,6) = Z-5(6,T)(/)(x,0) . 
Now we have additional renormalizations of and A 

ulib, T) = T)u^{b, T) Xnib, T) = Z\b, T)X{b, T) , 
and then it follows that 

du\ _ d{\og ^2 ^ ^-1 ^ 



db db \db 

dXn _ , rf(log Z) ^ ^^-JdX\ 



z=i 
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However, we have noted that to the order in which we are working, (1-loop), 
Z = 1. Since we shall be expanding in e this means that Z = 1 to 0{e) and 
so d(log Z)/db ~ 0{e^). Henceforth, we set Z = 1 above, and identify {uji,Xr) 

with (m, a). 

It is important to remark at this stage that the renormalization choices I am 
making (e.g. defining the coefficient of the kinetic term to be unity), are not 
forced on me by physics: they are convenient for the perturbative-style analysis. 
Other choices may be necessary in more complex situations in order to reveal the 
structure of more complex phase transitions in the most effective way. 
The transformation takes the form 

rr/j 2 A \ integration r-z^-i^ 2n\ n\ -6 a \ 

field renorm. r^r, 2n\ n\ -b \ \ 
— > S{(j),m {b),gn{b),e Aq) . 

Here Z = Z{b) and gn is the coupling associated with field monomial of order 
n. This forms the major part of the renormalization group transformation to be 
discussed later. The derived equations are a form of the renormalization group 
fiow equations for (m^, A). Note that by dealing with dimensionless quantities 
all explicit reference to the cut-off A disappears. 

Consider a fiow equation of the form 

This equation has fixed points A* where /3(A*) = 0, e.g.. 




If \*o < A(0, T) < A* then limj^oo A(6, T) = A* 

If A* < A(0, T) < A^ then limh^^^ X{b, T) = A* 

Fixed points A* where /5'(Ap) < are Infra- Red attractive, and alternatively if 
/3'(A*) > the points are Ultra- Violet attractive. 
Note that 

(1) In statistical mechanics we are given the UV or bare coupling A(0,T), 
and the IR or renormalized coupling is an IR point of /3(A). 
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(2) In quantum field theory we are given the IR or renormalized coupling 
A(oo,T): it is a renormalization condition from experiment or other low 
energy condition. Since A(oo,T) is not generally a fixed point of /3(A) it 
must be that the UV or bare coupling is a UV fixed point of (3. (It should 
be noted that although we have worked exclusively in Euclidean space it is 
believed that the divergences are the same as for the Minkowski version of 
the theory. This belief can be demonstrated in perturbation theory.) 

In our calculation we find the fixed points (w*-^, A*) 

e < (0,0) trivial f.p. - IR attractive 

e > (0,0) trivial f.p. - UV attractive 

( — , e) non-trivial f.p. - IR attractive 

It is believed that there are no others, not even at A = cxd. This means for 
quantum field theory that there is no UV fixed point to be associated with the 
bare coupling except the one at (0, 0) when e < 0. Hence, it follow that 

e < Xr =0 

e > Xr < if^e, 

where A^ = A(oo,T). This is the statement of triviality for quantum field 
theory. It means that in dimension four or greater the only consistent renormal- 
ized (p^ theory in the limit of infinite cut-off is a free theory: Xr = 0. Of course, 
we are not required to take the cut-off to infinity, rather we should keep it much 
larger than the physical momentum scales on which the theory is being applied. 
Even then, for finite A the analysis presented gives an upper limit to Xr. Such 
upper limits have been used to set upper bounds on the Higgs mass in the most 
common version of the standard model. 

Since e is small we can get a reliable expansion in e for both u*'^ and A*. This is 
the epsilon expansion 

oo 

A* = ^a„e". 

n=l 

It is hoped that the radius of convergence is greater than one, thus including 

results for D = 3. 

Write the equation for as 

We now show how to derive the postulated scaling behaviour from the fiow equa- 
tions. 

For «C 1 all parameters have reached their low-energy renormalized values. 
Let 

^ = > 1 b = log(^^Ao) . 
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We have 

This follows since for ^ ^ 1, m{b,T) = Then 



b-b 



lu^b,T) f{v?) 

= log(^eA) 

with A = e"*Ao and «C A «C Aq, but such that 

A(6,T) f« A* . 

This will always be possible if ^ is big enough. The key point is that we consider 
the flow first from Aq to A so that A flows very close to its IR attractive fixed 
point. We still have the condition that k ^ ^. We then consider the flow to 
A <C i^^. For this part of the flow A is fixed at its fixed point value A*. It is from 
tliis latter part of the flow that the scaling form for ^ as a function of t follows. 
The sequence of flows can be summarized thus 

Ao > A > >A 

A(0, T) A(6, T) f« A* A(6, T) f« A* 

0^6 ^ b 

Mow, for T = Tc, ^ = oo. Hence it follow from above that 

dv:" 



I' 



OO 



Thus u{b,Tc) = u* and the integral has a logarithmic singularity. We write 

T-T. 



u\b,%) = u*^ + Bt 



Then 



/ 



= log m + c 



di{t) Bdt 



at) f{u*^ + Bt) 

1 dt 



f'{u*^) t ' 



From above we have that 
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Hence we obtain the prediction that 

e(t) = Dt " with v = - + —e. 
^ 2 12 

The solution of the full flow equations looks like 



cutoff = A 



(u*^ g*) 




line of theories 
parametrized by T 



A, 



' 



critical surface 



/ 




intersection for T = 



We remark the following points 

(1) In a general description, the theory for some intermediate cut-off will be 
given by the values of a complete set of coupling constants associated with 
all possible field monomials in the action. These couplings, denoted hi, will 
be functions of A and T, and they flow in such a way that the low-energy 
predictions of the theory are independent of A. We would expect the flow 
equations to take the form 

The diagram of the flow is then, in principle, an oo-dimensional space. 

(2) The critical surface is the manifold of all theories in coupling constant 
space, h, which flow to a fixed point. If there is more than one fixed 
point then the critical surface is divided into domains of attraction each 
associated with one fixed point. All theories in the critical surface have 
^ = CO, i.e., they are at a continuous transition. This follows because it 
requires 6 — >• oo to flow to = which is the reference theory. Hence for 
all finite b we must have A = 6Ao > which implies that ^ = oo. 
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(3) As T changes the theory follows a trajectory (dotted line) which intersects 
the critical surface at T = Tc. As shown above the critical exponents are 
given by the behaviour of trajectories in the neighbourhood of the IR 
fixed point. Hence we only need to calculate /3j(h) in this neighbourhood. 

(4) The number of parameters we need to tune in order to intercept the critical 
surface is clearly equal to the number of unstable directions at the fixed 
point. In our calculation the number of such unstable parameters is one, 
namely the temperature T (i.e., u^{b,T) which is controlled by T). Thus 
the transition corresponds to an ordinary critical point. 

A fixed point with two unstable directions requires two external variables to 
be tuned e.g., T, P, and thus corresponds to a tricritical or bicritical point. 

(5) The critical surface separates the two phases i.e., the symmetric phase from 
the broken phase. 

(6) We expand A(h) about the fixed point, h*. This linearizes the fiow equa- 
tions. Writing A, = /ij - h* 

^ = R,,A, + O(A^) . 

The eigenvectors, e^, of Rij define the fiow directions and the associated 
eigenvalues, /x^ determine the stability properties 

fip < stable, /Xj, > unstable, fip = tricky . 



(7) For operator gn'P^"' we have [gn] = dn = n{D — 2) — D. Then /3rj(h) has the 
linear term 



^„(h) = -d^hn + . . . 

(The hn are the dimensionless versions of the g^.) This coupling will be 
stable for d„ > unless there are any very large renormalizations from 
other terms. Although large corrections cannot be ruled out it is expected 
that since, in practical cases, the dn are integers we find that the operators 
are stable if 

D 



Such operators are said to be irrelevant since their couplings automatically 
fiow onto the critical surface and do not need to be "tuned" . 

D = 4 n>2 . . . irrelevant 

D = 3 rj > 3 (p^, ■ ■ ■ irrelevant 



This refiects the result of our earlier perturbative calculation of Dc and 
the criteria for the validity of Landau theory. 
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(7) The flow equations are a coupled non-linear set. The hn of difl^erent fleld 
monomials get mixed. This reflects the operator mixing induced by renor- 
malization. In the linearized version the eigenvalues of Rij deflne eigenop- 
erators: linear combinations of field monomials which do not mix with each 
other. This can be generalized to the full non-linear version of the flow 
equations. 

(8) In the Landau theory we associate the coefficient of with the low- 
momentum dimensionful coupling 

^(T) = A^ lim e-*A(6,T) . 

0— >oo 

This is because g{T) is the coefficient of in the efFective action derived 
by the Legendre transform. = 1 to 0(e) otherwise we would need to 
include in the definition of g{T)). 

Let 

g{b,T) = e-''X{b,T) . 

Then the equation for the flow of A in the neighbourhood of the flxed point 
implies 



dg ^ 3_ 

db 167r2 



or 



and thus 



J L = ^ e^ft ( e{b--b) _ ^\ 



g{b) = 



g{b) 



1 + <f (e<*-^) - l) ■ 



Denoting g(T) = g{oo,T) we have the following results 

e > : g{T) = 0. Hence Landau's method fails and the relevant physics is 
encoded in how g{b, T) — >• as 6 — >• oo. 

e < : We have 

Since here g{T) is flnite and non-zero Landau's method works. As 
before we are assuming that b is large enough so that the equations 
we integrated to get this result are valid. In this case g{T) is actually 
independent of b. 
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: g{T) = 0. Landau theory fails but 

9{b) 



9{b) 



^^^^ ^ T+^lo.A- 
I.e., only logarithmic violations of Landau theory. 



4.5 The Renormalization Group in General 

Our procedure so far consists of "thinning" the degrees of freedom by integrating 
out the high momentum modes. We can define the transformation (A < 1) 

S{\)H{(j)) = H{\<l),\k) , 

with H{(j)) = H{l,<f),A) as the bare or microscopic Hamiltonian. The expUcit 
A dependence is there because the couplings have changed in a A-dependent way. 
The 5(A) form a semi-group 

S{M)S{\2) = S{M\2) . 

There is no inverse because of the existence of fixed points: information lost 
cannot be retrieved. 

It is convenient to perform a rescaling at this stage 

X — )■ x' = Ax length rescaling 

q — )■ q' = A^^q momentum rescaling 

<f> ^ (j)' <P'{fl') = ^^('^)^^(q) field rescaling 

AA — >• A cut-off rescaling 

The cut-off is rescaled back to its original value. The new fields (j)' have support 
< q' < A in momentum space, i.e., the same support as the original fields 
before thinning. A{X) is chosen by convenience. It may be chosen so that the 
coefficient of |(V^)^ in the new Hamiltonian is one, but other choices can turn 
out to be more appropriate in more complex theories. 
We then define the renormalized Hamiltonian Hji{<f), A) by 

nR{<l)',A) = H{\,<l),\A) . 

For example consider the kinetic term in H 

Change variables to q' = A^q 
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If we choose A{X) — Z{X)X (^+2), then the kinetic term in Hr{(I),A) has the 
canonical form. 

The renormalization group transformation is R{X) where 

R{X)H{(t>,A) = Hii{(t>,A) . 

Hr is defined on the same lattice as is H, but all dimensionful observables are 

scaled in the appropriate manner. 

Define 

G(qi, . . . , q„) =< (/)(qi) . . . (/)(q„) > . 

Then we must have 

< (/)(qi) . . . ^(q„) > = ^t(A) < ,^(31) . . . <^(^) >^ 
computed with H computed with Hr 

G(qi,...,qJ = ^?(A)G^(ai,...,^) 

With the choice for ^(A) above, Hr is a function of the set of dimensionless 
couplings which we denote here by g(A). then 

iJ^((/),A) = i?(g(A),(/),A) , 

and 

R{X)H{so, A) = i?(g(A), (/), A) , go = g(l) . 

This defines the RG transformation on the couplings 

RW go = g(A) 
and since R{Xi)R{X2) = R{XiX2) we have 

^(A2) g(Ai) = g(AiA2) . 

Thus 

R{1 - a) g(A) = g(A) - «A^g(A) + . . . 
A(g) = 

a log A 

= \im -[R{l-a)gi{X)-gi{X)] . 

This defines the function A(g)- 

(1) /3i has no explicit dependence on A if g is a vector of dimensionless couplings. 
The gi correspond to the dimensionful couplings scaled by A. 
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(2) /3j does not depend explicitly on A: the only A-dependence is through g(A) 
in H. This follows since R{1 — a) is independent of A. The next step along 
the flow in coupling constant space cannot depend on A since A tells us 
how far we have already evolved, or flowed, along the trajectory. However, 
the next step depends only on the couplings at that point on the flow and 
not on the flow history. This is guaranteed by the semi-group property of 
the RG transformation. 

Setting b = — log A we have 

dgi r,f ^ 

A fixed point is determined by A(g*) = 0, Vi Then 

i?(A)i?(g*,(/),A) = i?(g*,A) , 
and hence from above, when g = g* 

G(qi,...;q„) = (A)G(^, . . . ; ^) 
^(A) = X-^^'+^^ZiX) . 

Then writing Z{X) = A^''* we have 

^(A) = A-^(^-^^) , A, = - l) + , 

and thus 

G(qi,...;q„) = A-"(^-^^)G(^, . . . ; ^) 
This is the scaling law for Green functions at the fixed point. 



A(f, is the scaling dimension of (/)(x) 
— 1^ is the engineering dimension of (ji(x) 
r]^ is the anomalous dimension of ^(x). 

The full renormalization group thus consists of two stages. 
Momentum space 



A AA A 

q q A^q 

r escaj e 

m m X^^m 

^(q) ^(q) ^'(q') = ^-^^(q) 

= A(^-^*)(/.(q) 
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Coordinate space 



a a/A a 

X X Ax 

'S'(A) rescale 



(/>(x) (/.(x) (/.'(x') = A-^^-^<^(x) 

(1) The special form for <f)' applies at the fixed point. 

(2) Since the full transformation, R{X), leaves A, or a unchanged the overall 
effect is a rescaling of the physical parameters and observables only. 
A plays no role in the dimensional analysis based on this rescaling. 

(3) A general field combination /(^(x)) generally "renormalizes" in a compli- 
cated way 

m 4 /'((/)) , 

but for special combinations, the scaling fields, a dimension, Af, can be 
assigned 

/'(<^(x)) = A^//(<^'(x')) . 
This is possible only near a fixed point. 

(*) for information only - not directly part of course 

(i) Before thinning have term in action of the form 

/ d^'x /i/((/)(x)) , 

J a 

where /i is a coupling or external field. After thinning this term becomes 

/ d^'x hf'{<l){^)) . 

Ja/X 

Near to the fixed point we can write this term alternatively as 

/ d^'x hX^ff{4>'{^)) , x' = Ax . 

./a/A 

= / d^x'hX^f-^f{<i^'{^)). 

J a 

But from the RG we are told that this term in the rescaled action must be 

/ d^'x hx-^^f{4>{^)) . 

J a 

and hence we verify the existence of a scaling dimension for the operator 
function / with dimension Aj which, comparing the last two expressions, 
is given by A/ = D — A/j. 
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(ii) An alternate way of seeing this is to use the scahng rules for Green functions 
near to a fixed point. These rules are derived in section 5 and you should 
refer back to here after reading it. We have 

d f 

— G(xi, . . . , x„; /i) = j d^y G/(xi, . . . , x„, y; /i) , 

where Gj means an insertion of /(^) into the expectation value. By the 
association of a dimension Ay^ with / we use the scaling rules to write this 
in the alternative form: 

= j d^'y A"^*+^/ Gj(Axi, . . . , Ax„, Xy- X-'^'^h) . 
However, this last expression must also be given by 

|^[A"^^G(Axi,...,Ax„;A-^''/i)] 
= j d^'y X^'^^X-^^ G;(Axi, . . . , Ax„, y- X-^^h) 
= A^ y d^'y' A^^^A-^'' Gj(Axi, . . . , Ax„, Xy'- X'^^h) . 

Comparing these two versions we find Af = A — Ah as before. 

(4) The anomalous dimension, 77,^, is computed from Z{X), the field renormal- 
ization constant. 

Now consider the flow in the neighbourhood of a fixed point. Write 

* I 

g = g + V . 

Then 

dvi 



db 



The eigenvalues of Kij are with eigenvectors e^. 
Then we find 

a 

= Aa{g*)ha no sum on a, 
ho. = /^L^^A-^". 
Remember that b — log A , A = A/Aq < 1. 

If Aq > the eigencoupling ha is unstable and is repelled by the fixed point 
and vice- versa if A^ < 0. We have recovered the result that we stated before 
in a slightly different context namely that every fixed point corresponds to a 
continuous transition and can be labelled by the number r of repulsive (i.e.. 
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positive) eigenvalues, then r is the number of external fields that must be tuned 
so that we see the transition controlled by that fixed point. 

In general the Hamiltonian is 

H = fdxY^ 9nAn{x) = H*+ f dxY^ VnAn{x) , 

where the A^ are general operators and are functions of the fields. Near a fixed 
point we can write instead 

H = H*+ [d^Yl ^a^a(x) , 

where 

= T,n (ea)ri^ri(x) SCaliug fieWs 

ha = {Ua)n9n SCallug COUpliugS. 

Here we have defined the vectors Uq by 

Now we define [Aa] = A^^ and we have that [ha] = Aq. But H is dimensionless 
and so we must have 

AA^=D-Aa. 

Aa is an eigenvalue of the linearized RG transformation. 
Denote Green functions by 

G(xi, . . . , x„; h) =< ^(xi . . . ^(x„) >c , 

where h is the vector of eigencouplings. Near a fixed point we have 

G(xi, . . . , x„; h) = it(A)G(Axi, . . . , Ax„; A-^"h) , 

where A^ = X^A^. The extra factor of A^ occurs when G is fourier transformed 
from momentum space to coordinate space. But we have 

i3 = A^^3 = , 

and hence it follows that 

G(xi, . . . , x„; h) = (A)"^^G(Axi, . . . , Ax„; A-^"h) . 

This is the full scaling result for Green functions evaluated near to a fixed point. 
In the next section we shall derive the scaling theory which demonstrates that 
such results can be obtained even when the bare theory is not near to a fixed 
point as long as it is sufficiently close to a critical surface, i.e., near a continuous 
phase transition. 
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Note for example that the "magnetic field" , h, has dimension 

since the relevant term in the Hamiltonian is — J dx hcj) which must be dimen- 
sionless. 

If Aa is a relevant operator then we must have 

A,,^ = D - > ^ ^A^<D . 
So, for example, ^ is a relevant operator. 

We can see an immediate application of these ideas as follows. Let 5(x) have 
dimension which does not necessarily satisfy A^ < D, i.e., B need not be 
a relevant operator. Let ^ be a relevant operator conjugate to external field 
h: A^ = D — Ah and A^ < D. All other parameters except h are fixed at 
their critical values, and so as /i — )■ ^ diverges. The corresponding generalized 
susceptibility is 

^ = Jl[d^, <S(OM(xi)...^(xJ>, . 
[This follows since interesting term in i? is —/i/ fix ^(x) and we have 

< B{0) >= j del) e-^ B{(l){0)). 

] 

The only length scale in the problem is ^ and hence on dimensional grounds 

and for n large enough this will diverge as h ^ 0, ^ — )■ oo if 

n{D - Aa) > Ab . 

This will occur for some n since A is relevant and hence A^ < D. 
Thus singular thermodynamic quantities are associated with relevant operators 
and consequently with the unstable or repulsive eigenvalues of the linearized 
RG transformation. 

5 Scaling theory 

Consider the Green function for a field ^(x) which has dimension A<^ 

G'(xi, . . . , x„; g|^) =< (j){-Xi) . . . (j){-Xn) >c ■ 

G has been labelled by ^ which will be assumed to be large. ^ is a dynamical 
variable determined by g, the couplings in dimensionless form. From the RG we 
have, in general 

G(xi, . . . , x„; g|0 = i^(A)G(Axi, . . . , Ax„; g(A)|AO ■ 
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Now consider a particular value, A, of A such that for cases where ^ 1 

A < 1 and A^ » 1 and g(A) = g* . 

This will always be possible if ^ is large enough, although how large it must be 
for this to be true is a question we have not addressed. Then we can define the 
function F by 

G(xi, . . . , x„; g|0 = F{xi, . . . , x„|0 ■ 

Which, from above, means that 

F(xi, . . . , x„|0 = it(A)G(Axi, . . . , Ax„; g*|AO ■ 

Then a further renormalization gives 

F(xi,...,x„|0 = i5(A)F(Axi,...,Ax„|A0 
= A"^*F(Axi,...,Ax„|AO . 

The last result follows since this further renormalization is applied to a theory 
near the fixed point (i.e. only the mass, m = is not at its critical value) and 
so A takes the power-law form dictated by the linearized RG flow equations. 
Thus for = p, 1 and fixed we have 

G(xi,...,x„;g|0 = r"^*£^(6,---,U , 
where = Xj/^ and 

= p"^Mt(A)G(Axi,...,Ax„;g*|Ap) . 

This result is a verification of the scaling hypothesis. As can be seen, all 
results are ultimately expressed in terms of a theory with correlation length \p. 
In order that this theory lies in that neighbourhood of the fixed point in which 
the linearized RG equations apply we must have that Xp is sufficiently large. A 
diagram of the fiow is shown below. 



within linear 
region of RG 



correlation length line of theories 



I 



g 



What matters is that the theory with correlation length Xp can be so chosen to 
lie in the region where the linear RG approximation is good. This region is a 
neighbourhood of the fixed point and as long as the fiow line emanating from the 
bare theory passes through this region then the above analysis can be applied and 
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the scaling behaviour derived is vahd. Clearly, the range of values of temperature 
T (near Tc) for which this happens depends on how close the intersection of the 
trajectory of bare theories with the critical surface is to g*: if it is close to g* then 
the scaling behaviour will be apparent over quite a wide range of temperatures 
near Tc, but the range will narrow as the distance of the intersection from g* 
increases. This range of temperatures in which the scaling behaviour is observed 
and the critical exponents can be measured is called the critical region. 

Now follows a few interesting consequences of scaling theory. Consider the gen- 
eralized susceptibility 



r„ = J dxi... dxn-i G(0, xi, . . . , x„_i; g|0 



We have 



" dh dh d^ 

where h is conjugate to <f). 
From above 

where C„ is a constant given by 

Cn = J C?%...C?%_iE(0,6,...,^.-l) . 

Then from above 

r„ = [{n - 2)D -{n- l)A^]C„_i ^(-2)^-(-i)a,-i ^ 

dh ^ 

1 

^ ~ /i""'* where z^, 



Thus the critical exponent associated with the relevant external field h conjugate 
to (j), is given in terms of the scaling dimension of ^, namely A<^. 

5.1 Scaling Relations 

(1) From above the magnetic susceptibility satisfies the relation 

Since 7 is defined by x ~ t^'^ we have the relation 
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(2) It follows from scaling theory that 

M =< (/) >~ C^i' ~ . 
Since j3 is defined by M ~ t'^ we have that 

(3) /3 (= 1/T) is conjugate to the energy operator e(x) and hence 

^ ~ t^" with v 



D-A, 

Strictly speaking the only part of e(x) that matters is the most singular, 
i.e., the component that has the lowest scaling dimension. In theory this 
is the operator <p^. 

(4) The specific heat is 

'dS\ _ fd'^F\ 



V 



with F = -Tlog Z. Thus 

Cv = fpo,Z 

= Jd^ < e(0)e(x) > ^ 

Since a in defined by Cy ~ i^" we find 

(5) From above M ~ ^^'^"^ and ^ ~ /i . Since S is defined by M ~ /i« we 
find 



S = 



A^ 



Combining these relations we derive the scaling relations 

a + 2^ + 7 = 2 
a + 2(35 = 2 + 7 . 

Note the following 

(i) Only two relevant or unstable directions occur in 0^ field theory which 
correspond to h and T and are associated with the relevant operators (j) and 
e respectively. 

(ii) All critical indices are controlled by A<^ and or equivalently by the di- 
mensions assigned to h and T 

[h] = D- A^ and [T] = D - A, . 

The rest is dimensional analysis. Hence in this case there are only two 
independent indices. 
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5.2 Renormalized Green Functions 

This is a very large topic since the ideas underlying the concept of renormalization 
and their implementation are central to the modern theory of analyzing both 
quantum field theory and quantum and classical statistical mechanics. In this 
course we have time only to introduce the main ideas using the critical theory as 
the example, and to link these ideas with the notions which have gone before. 
We work with a scalar field theories which lie in the critical surface 
= oo = 0). Of course, the approach can be extended to study 

near critical theories and non-leading corrections to the predicted leading scaling 
behaviour. 

Define r„(p, go. A) to be the IPI truncated Green function. We will consider the 
effect of the "thinning" procedure but will not rescale the observable quantities. 
In other words we shall use the fact that the "thinning" procedure produces a 
new theory with the same low-energy predictions as the original one. The effect 
of "thinning" is encoded directly in the relation 

r„(p,go,A) = Z-i l^^j r„(p,g(/x),/x) . 

/i is the new cut-off after thinning and the g(/x) are dimensionless couplings with 
g(A) = go- 

The LHS is independent of /i and hence 

(^//|^ + A(g)^ - 7^(g)) r„(p,g(/.),/.) = , 

where 

d 

A(g) = l^-Q^3i{lj) ■ 

This equation satisfied by r„ is the renormalization group equation. 

Note that 

lim r„(p,g(/x),/x) = lim ( - ] r„(p,go,A) 

A— >oo A— >oo \ 1^ / 

is finite. r„(p, g(/i), is the renormalized Green function. 
Consider field theory. 



(i) 7(5) 



^r2(p,^o,A)p2=^2 = z M ^r2(p,5(/x),/x)p2=^2 



= (^) KM) . 
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f{g) is a dimensionless function of g and can be calculated in perturbation 
theory 

f{g) = l + ag^ + ... 

There are no divergences because |p| = /i, which is the cut-off. We also 
have (calculating in the limit e — >■ 0) 

^r2(p,^o,A)p2=/x2 = l + 5olog^^j +■■■ 

The corrections introduced by keeping e finite in this calculation give rise 
to corrections to 7 which are higher order in e. Then from above we get 

^/A\ ^ l + ag\fi) + ... 



./^y l + agi\og[^) + ... 
d 

^ 7 = log Z 

O/J, 

= ago +2ag{n)(3{g{n))... 

From the renormalization group equation 7 is a finite (i.e., A independent) 
function of 5'(/i). This follows because both I3{g) and Tn{p, g{iJ.), i-i) are in- 
dependent of A and the RG equation cannot balance unless 7 is independent 
of A too. 

Since ^'o = fl'(/^) + 0{g^{ii)) we can make the substitution for g^ in terms 
of g{jjL) and find 

7(5) = ag" + 2agM ■ ■ ■ 
Ti{p,go,A)p2=^2 = Z'"^ r4(p,^(/x),/x)p2=^2 . 

But 

r4(p, g{^J■), /^)p2=,i2 = ^J■\g{^J■) + bg^{^i) + ...) , 

again a finite function: b is an e-independent number. Note that r4 has 
engineering or ordinary dimension e: [r4] = e. These dimensions are 
carried by the /x^ factor. 

We also have that 

r4(p, go, A)p2=;.2 = A^^o - ^A2>-^^2 _^ 

To see how the form of the second term on the RHS arises we note that it 
comes from an integral of the form 



'■7 



Const. , 
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Ignore the /i-independent terms and terms that are not leading non-trivial 
order in e. 

[We note in passing that the above equations show to lowest order that 
g{^) = (^j^ go- In the calculation for 7 we worked in the e — )■ limit since 

we only required the leading order results in that case. This is why we were 
able to use go — g{fi) in that calculation. If we want to work to higher 
orders in e we must be more careful.] 

Then 

+ bg^fi)) ^ A^go - ^A'^-^^o , 

(Note that Z = 1 to the order in which we are working.) Now apply the 
operator ji^ with /3 = n^^. We get 

+ + W + 2bg(3) f« (3o(^^^ gl 

= /3o/ + ..., 

where the expression for go in terms of g has been used in the last step. 
Hence 

_ -eg + {/3o - eb)g^ 
^~ 1 + 2bg 

Keeping only non-trivial orders in e, (/3o is gives 

^ = -eg + /3o/ + . . . 

We have already computed /3o before when calculating the critical exponent 
u (note that when comparing these two calculations we must remember 
that (ilog/i = —db. This accounts for the sign difference between the two 
otherwise identical versions of the flow equation). For /i -C A we find 
g{li) = g* where 

W) = ^ 9* = ^. 

167r2 . 



and we found before that g* = ^^^e. 



A consequence of the two calculations above is that for «C A we have 

^ z(-]=j{g{fi)) = j{g*). 



dlog/i \fj, 
Or, writing A = ^ 
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But from the scaling analysis near the fixed point 
Thus 



V<t> = 








= 




since 


Kg*) = 



We now indicate how to calculate a. T2 is given by the graphs 




)"' + p + p >f 3- + 




p 

+ CA'gl + (5oW + ap^)log- +... 

The calculation of a is somewhat involved and will not be reproduced here. From 
Raymond's book "A Modern Primer on Field Theory" we find 

1 

a 



6(167r2) 



2 



1 1 (]^\ 2 

'^'^ ~ 2 6(167r2)2 1^~3~ J ^ 

= 

108 



What happened to all the couplings corresponding to the irrelevant operators 
(i.e., those with dimension A > D)? A calculation of the /3 functions for these 
operators in the same way as in (ii) above shows that the fixed point for these 
couplings is at the origin. This is expected because 

(i) the operators corresponding to these couplings do not give rise to extra IR 
divergences: D > Dc for these operators. This implies that Landau's ap- 
proach works and the critical exponents receive no anomalous contributions 
from these operators. 

(ii) mean field theory works for operators whose associated couplings have their 
fixed point at the origin. 

Vci> = ^7(g*) , 
but only g*, the couphng, is non-zero. 
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(iii) The corollary is that only the renormalizable operators in the usual sense, 
i.e. A < D are relevant since only their couplings flow to IR stable 
fixed points which are not at the origin. For those operators with A = D 
the fixed point is at the origin but we find log, rather than power-law, 
corrections to mean field theory. These operators constitute a special case 
but an important one since, for example, scaling violations in Quantum 
Chromodynamics are exactly of this kind. 



